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Abstract: We study dynamical aspects of the plane- wave matrix model at finite temper- 
ature. One-loop calculation around general classical vacua is performed using the back- 
ground field method, and the integration over the gauge field moduli is carried out both 
analytically and numerically. In addition to the trivial vacuum, which corresponds to a 
single M5-brane at zero temperature, we consider general static fuzzy-sphere type config- 
urations. They are all 1/2 BPS, and hence degenerate at zero temperature due to super- 
symmetry. This degeneracy is resolved, however, at finite temperature, and we identify 
the configuration that gives the smallest free energy at each temperature. The Hagedorn 
transition in each vacuum is studied by using the eigenvalue density method for the gauge 
field moduli, and the free energy as well as the Polyakov line is obtained analytically near 
the critical point. This reveals the existence of fuzzy sphere phases, which may correspond 
to the plasma-ball phases in A/" = 4 SU (oo) SYM on x S^. We also perform Monte 
Carlo simulation to integrate over the gauge field moduli. While this confirms the validity 
of the analytic results near the critical point, it also shows that the trivial vacuum gives 
the smallest free energy throughout the high temperature regime. 
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1. Introduction 

Matrix models provide a promising approach to non-pcrturbative dynamics of superstring 
theories and M-theory [1-3]. In particular the matrix model proposed by Banks, Fischler, 
Shenker and Susskind (BFSS) [1] is conjectured to be a non-perturbative formulation of 
M-theory, a hypothetical theory in eleven dimensions, whose low energy effective theory is 
given by lid supergravity. The BFSS matrix model, which takes the form of a matrix quan- 
tum mechanics, can be obtained also from the supermembrane theory in eleven dimensions 
via the matrix regularization [4]. Studying the BFSS matrix model at finite temperature is 
an interesting subject due to its relation to the Schwarzschild black hole [5-7]. Such studies 
have been performed, for instance, in refs. [8,9]. Applying matrix models to cosmology is 
another fascinating direction, which has developed recently [10]. 

While the BFSS matrix model is defined on a flat Minkowski space-time, we may 
consider matrix models also on curved space-time. In this regard the so-called pp-wave 
background provides a particularly simple and tractable example. In fact the maximally su- 
persymmetric pp-wave background in eleven dimensions, which preserves 32 supercharges. 
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is unique and its explicit form is known [11]. Berenstein, Maldacena and Nastase [12] pro- 
posed a matrix model on this background, which is closely related to the supermembrane 
theory on the pp-wave through the matrix regularization [13-15]. It has fuzzy spheres as 
classical solutions due to the presence of the mass term and the Myers term [16]. The 
stability of these solutions [13, 17, 18], as well as the spectrum around them [13, 19-21], 
has been studied intensively, while more general classical solutions such as a rotating fuzzy 
sphere are also discussed [18,22-26] ^. In fact the fuzzy sphere solutions can be interpreted 
as giant gravitons, and the interaction potential between them ^ is shown to be of the 
1/r^-type [18,26] similarly to the results in the BFSS matrix model. This type of po- 
tential is anticipated since the spectrum in the linearized eleven-dimensional supergravity 
around the pp-wavc background [32] is included in the matrix model as the spectrum of 
the zero- mode Hamiltonian [15,20]. 

The plane-wave matrix model has been studied also at finite temperature by various 
authors. In ref. [33] the free energy around the trivial vacuum, which corresponds to a 
transverse M5-branc [34] at zero temperature ^, was evaluated at the one- loop level, and 
the Hagedorn transition was studied in detail. (See refs. [36, 37] for a two-loop extension 
and ref. [38] for a review on this subject.) 

In this paper we extend such calculations to general classical vacua, and identify the 
vacuum that gives the smallest free energy at each temperature. Similar studies have been 
made earlier in refs. [39,40]. However, the degrees of freedom corresponding to the gauge 
field moduli, which appear at finite temperature, were not taken into account although they 
actually play an important role in the Hagedorn transition [33] . Here we perform the one- 
loop calculation by the background field method keeping the gauge field moduli arbitrary. 
The integration over the gauge field moduli is performed both analytically and numerically. 
Near the critical point we use the eigenvalue density method to obtain analytical results 
extending the previous works [33, 38] for the trivial vacuum. This reveals the existence of 
fuzzy-sphere phases, which may correspond to the plasma-ball phases in A/" = 4 SU(oo) 
SYM on X S^. We also perform Monte Carlo simulation to integrate over the gauge 
field moduli with the matrix size up to = 2000, which allows a reliable extrapolation 
to = oo. (Note that the results of ref. [39,40] are obtained for small N.) While this 
confirms the validity of our analytical results, it also shows that the trivial vacuum gives 
the smallest free energy throughout the high temperature regime. In fact the leading 
asymptotic behavior in the high temperature limit is universal, but the difference appears 
in the subleading term. 

Comparison of the free energy for general classical vacua was also made in a matrix 
model, which is obtained by dimensionally reducing the plane-wave matrix model to a 
point [41]. There the true vacuum turned out to be a complicated multi-fuzzy-sphere type 

matrix string theory on a type II A plane- wave background, which includes fuzzy spheres as classical 
solutions, has originally been constructed in [27]. The spectrum around the fuzzy spheres was computed 
in [28] by following the method of [13]. The theory is applicable to the matrix cosmology scenario [29]. 

^The interaction potential between point-like gravitons has been studied recently in ref. [30]. See also 
ref. [31] for a short review on the graviton potential. 

^In a deformed plane- wave matrix model with an interaction term coming from the 6-form potential, a 
fuzzy five-sphere solution was constructed [35]. 
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configuration, wliicli incorporates a non-trivial gauge group. Analogous models without 
the mass term have been studied intensively by Monte Carlo simulation [42] . These studies 
are partly motivated from the IIB matrix model [2], a conjectured nonperturbative for- 
mulation of type IIB superstring theory, which can be formally obtained by dimensionally 
reducing the BFSS matrix model to a point. More intimate relationship between the two 
models has been suggested [43] on account of the large N reduction [44]. Calculation of the 
free energy in the IIB matrix model for various configurations provided certain evidences 
that four-dimensional space-time is favored dynamically [45] . See ref . [46] for related works 
on this issue. 

The results obtained in the present paper are expected to have close relationship to 
the phase structure of large N gauge theories at finite temperature and at finite volume, 
which has been explored in refs. [47-51]. According to the AdS/CFT correspondence [52] 
at finite temperature [53], the observed phase transitions should correspond to those on 
the gravity side such as the Hawking-Page phase transition [54] and the Gregory-Laflamme 
phase transition [55]. This provides an explicit realization of the connection between the 
Hagedorn transition in string theory and the deconfinement phase transition in large N 
gauge theories pointed out earlier [56]. 

The rest of this paper is organized as follows. In section 2 we briefly review the plane- 
wave matrix model and its classical vacua. In section 3 we perform one-loop calculation 
around general classical vacua keeping the gauge field moduli arbitrary. In section 4 we 
perform the integration over the gauge field moduli to obtain explicit results for the free 
energy and the Polyakov line. Section 5 is devoted to a summary and discussions. In the 
Appendices we present the details of our calculations. 



2. The plane-wave matrix model and its classical solutions 

In this section we define the plane- wave matrix model and discuss its classical solutions. 
In order to study the model at finite temperature, we make the Wick rotation t — > —it and 
compactify the imaginary time direction to a circle with the circumference (5 = 1/T. Thus 
the action can be written as ^ 



/3 

5pp = I dtTi 





2 V3 



{X,f + - (^) (Xa)^ + if^eijkXiXjXk + Z^*t^l23* 



(2.1) 



where Df = dt — i[A, ■ ] represents the covariant derivative, and the indices M,N run 
from 1 to 9. Since the transverse 5*0(9) symmetry is broken down to 5*0(3) x 50(6) on 
the pp-wave background, we have also introduced the S*0(3) indices k = 1,2,3 and the 



*Note that the model studied there is closely related, but not equivalent to the high temperature limit 
of the plane- wave matrix model. 

®We have made the rescaling, A — » RA, t j^t , — » -R/u, in the action presented, for instance, in 
ref. [14] to arrive at the present form, which does not include the parameter R (the compactification radius 
of the 11th direction in M-theory) explicitly. 
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so (6) indices a = 4, • • • , 9 . The partition function for the finite temperature system is 
defined by 

Z = J[dAit)][dXMmd^md^Ht)]^^Pi-Spp) , (2.2) 

where the bosonic and fermionic fields obey the periodic and anti-periodic boundary con- 
ditions, respectively, in the imaginary time direction. 

Let us discuss static classical solutions of the model within the usual Ansatz Xa = 
^' = . At zero temperature it is convenient to set the one-dimensional gauge field A{t) to 
zero by choosing the gauge, while at finite temperature we can only set A{t) to a constant 
diagonal matrix due to the nontrivial holonomy in the temporal direction. In the case of 
the trivial vacuum Xj = (i = 1,2,3), for instance, the gauge field can be an arbitrary 
constant diagonal matrix without increasing the action, and the corresponding moduli 
parameters should be integrated. These degrees of freedom, which we refer to as the gauge 
field moduli, indeed play a crucial role in the Hagedorn transition [33] . 

In fact the plane-wave matrix model has static fuzzy-sphere type classical solutions 

X, = 5,^|l, = |0(l1-)01,,), (2.3) 

s 

Ait) = A{t)^^(ln,^A('Ht)) , (2.4) 
1=1 

where L^"^ are the SU{2) generators in the n-dimensional irreducible representation satis- 
fying = ieijkL^^^ , and the parameters kj and n/ obey nj ■ kj = N . The 
kj X ki matrices A^^\t) represent the gauge field moduli. One can easily check that the clas- 
sical action for the above solution is zero for arbitrary A^^^t) by noting that [i?*, ^(t)] = . 

At zero temperature, these solutions arc interpreted as stacks of M2-branes or M5- 
branes (or their mixture) depending on how one takes the large N limit. In what follows 
we often restrict ourselves to the s = 1 case for simplicity, and use the parameters n = n\ 
and k = ki, which satisfy n ■ k = N. Then the solution represents k copies of spherical 
M2-brancs if k is held fixed, while it represents n copies of spherical M5-branes if n is held 
fixed, in the large N limit. An evidence for the latter interpretation is provided by the 
coincidence of the spectrum [34] for the trivial vacuum (n = 1) case. Note, however, that 
the argument for identifying transverse M5-branes in the plane- wave matrix model is based 
on the protected BPS multiplet, and it relies rather crucially on supersymmetry, which is 
broken at finite temperature. Hence it is not totally clear whether the M5-branes are 
stable against thermal fluctuations. One way to demonstrate the existence of M5-branes is 
to compute their radius at finite temperature by the Gaussian expansion method extending 
the calculation in ref. [7] to the present case. 

3. One-loop effective action 

In this section we compute the one-loop effective action around the classical solutions 
keeping the gauge field moduli arbitrary using the background field method. The matrices 
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are decomposed as 



Xi{t) =Bi + Yi{t) , Xa{t) = + Ya{t) , 

A{t) = A{t) + A{t) , ^{t) = + -^{t) , (3.1) 

where Yi{t), Ya{t), A{t) and ^'(i) represent the fluctuations around the background Bi and 
A{t). 

For the trivial vacuum Xi = , since the fluctuation A{t) can be totally absorbed into 
the gauge field moduli A{t), we set A{t) to zero. Similarly, in the general case, we omit 
S/=i(^/)^ zero modes in the fluctuation A{t), which corresponds to changing the gauge 
field moduli. In the presence of the background, the U {N) gauge symmetry is broken 
down to 11/=! U{kj). Therefore, we also obtain zero modes in the direction of the gauge 
orbit corresponding to the broken symmetry. These zero modes should be removed by an 
appropriate gauge fixing, and we use the background field gauge for this purpose. By using 
the unbroken gauge symmetry, the gauge field moduli A^^^ (t) can be brought into a static 
diagonal form 

^«(0 = -^diag(af\---,4^;), (3.2) 

where € ^r] . See Appendix A for the details of the two gauge fixing procedures. 

In the rest of this section, we integrate out the fluctuations perturbatively keeping 
the moduli parameters (o = 1, ■ ' ' i arbitrary. We add the gauge flxing term S'g.f, 
and the ghost term Sghost given by eq. (A.6) to the action, and expand it with respect to 

the fluctuations as S = S'^o) + S^^^ H h S'(^). The classical action S^^^ vanishes for the 

backgrounds we consider here, and the linear term S^^^ vanishes, too, since the backgrounds 
satisfy the classical equations of motion. The quadratic term is given by 



= dtTr 
Jo 



^{DtYiY - [y,,y,-][B,,5,] - -[K„S,f + -[^) Yi 
-if^eijkYi [B,,Yj] + ^(DtYaf - ^[Fa,^^' + ^ 

+^(DtAf - Bi]^ + Dtc ■ Dtc - [B^, c][Bi, c] 



2 



(3.3) 



where c and c are the Faddeev-Popov ghosts. We have decomposed the spinor indices 
according to SU{2) x SU{4), where a = 1, 2 and A = l,--- ,4 represent SU{2) and SU{4) 
indices, respectively [13]. The contributions of the higher order terms can be neglected in 
the ji ^ oo limit, as can be seen by rescaling of the variables as 

Ym f^-^^^YM , A ^ , A^fiA, 

c fi-^l'^c , c ^r^l'^c , t [TH , (3.4) 

which brings the action S into the form 

5 = 5(2)+/,-3/25(3)+^-3^(4) ^ 
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where /i that appear in S^'^^ , S^^^ and S^^^ are set to unity. Hence the one-loop calculation 
is justified in the large ji limit. Since we also take the large A'^ limit, it is important to 
clarify the N dependence of the expansion parameter. Let us recall that the expansion 
parameter in perturbation theory is given by = (^^-i-yj. , where we have temporarily 
restored the parameter R. Clearly the expansion parameter diverges in the large N limit 
with fixed n and p'^ = N/R. Hence, in order to ensure the validity of the perturbativc 
expansion, we need to take the — > oo limit faster than the large N limit. Note also 
that non-perturbative effects such as the tunneling between the trivial vacuum and the 
fuzzy sphere vacua is suppressed due to the iV — oo and /x — oo limits [57]. With this 
understanding, we set n = 1 to simplify the expressions in what follows, but one can restore 
the fi dependence of the results after integrating the fluctuations by simply replacing /? by 
Pfi (i.e., replacing T by ^). 

Let us introduce the operators 

CiM=[Li,M], (3.5) 

which act on an N x N matrix M. Then the quadratic term S^"^^ in the action can be 
rewritten as 



r0 

5(2) = dtTr 
Jo 



^{DtYif + ^YiPi,{C)Y, + \{DtYaf + \YaQ{C)Ya (3.6) 



+V'^^'^AV'Aa + 4^^^''R^ai^)4^Aa + ^(A^)' + ^ir(£)i + DtcDtC + cT{C)c 

where we have defined the following mass operators 

Pij{C) = ^ { [Ck'' + 1) 5ij - ieijkCk} , Q{C) = ^ [Ck"" + , 

R^ai^) = \€ + \4a^k , T{C) = ^£,2 . (3.7) 

The mass spectra, given by the eigenvalues of these operators, have been obtained for 
the physical modes [13] and for the unphysical modes [18]. Due to the structure (2.3) 
of Lj, it suffices to solve the eigenvalue problem for the (n/ x n/) square block and the 
(n/ X nj) rectangular block in the matrix, on which the mass operators act. The results 
are summarized in tables 1 and 2, respectively. Note that the mass spectra for the square 
block in the gauge and ghost fluctuations do not contain massless modes unlike the spectra 
obtained in rcf. [18] for the vanishing background gauge field. In fact those massless modes 
are treated exactly as the gauge field moduli and an appropriate gauge fixing in the present 
formulation. 

Integrating out the fluctuations, we obtain the determinant of operators appearing in 
the quadratic action (3.6), which can be evaluated in a standard way, once the mass spectra 
are given. (See Appendix B.) Including the Vandermonde determinant, which comes from 
the gauge fixing for the moduli integration (See Appendix A. 2.), the partition function is 
obtained at the one-loop level as 

Z = C J [da] exp{-Ses[a]) , (3.8) 
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type of fluctuations 


mass 


spins 


degeneracy 


Yi (i=l,2,3) 




l<l<ni -1 


21 + 1 




< I <ni - 2 


2Z + 1 




¥ 


l<l<ni 


2Z + 1 


ya(a = 4,... ,9) 




0<l<ni-l 


6(2Z + 1) 


il) (fermion) 




l<l<ni-l 


4(2Z + 1) 




ii+T-2) 


i < Z < n, - i 


4(2Z + 1) 


A (gauge) 




1 < Z < n/ - 1 


2Z + 1 


c, c (ghost) 


3\/K« + l) 


l<l<ni-l 


2Z + 1 



Table 1: Mass spectrum for the (nj x m) square block. 



type of fluctuations 


mass 


spins 


degeneracy 


Yi(i= 1,2,3) 




i n/ - nj\ < I < i(n/ + nj) - 1 
i|ni - njj - 1 < / < i(ni + nj) - 2 
i|ni — nj + 1 < Z < i(n/ + rij) 


2Z + 1 
2Z + 1 
2Z + 1 


Fa(a = 4,--- ,9) 




i nj — nj\ < I < i(ni- +nj) — 1 


6(2Z + 1) 


ip (fermion) 


(i + i) 
(1 + 1^) 


^ - nj \ - 1 < Z < |(ni + rij) - | 
i n/ - nj + i < Z < i(n/ + nj) - i 


4(2Z + 1) 
4(2Z + 1) 


A (gauge) 




I n/ - nj < Z < i(n/ + nj) - 1 


2Z + 1 


c, c (ghost) 




i|n/ - nj| < Z < i(n/ + nj) - 1 


2Z + 1 



Table 2: Mass spectrum for the {m x nj) rectangular block. 



where -Sefffa] represents the effective action for the moduh parameters aa \ and the inte- 
gration measure [da] is given by 



1 t4- daa 



^ ^ ki\ J-J- AttT 

1=1 a=l 



The normalization constant C turns out to be the same for all the classical solutions, as we 
demonstrate in the Appendix C, and therefore it is irrelevant. 

For the trivial vacuum, the masses are given by ^ , g and | for the fluctuations Yi, Ya 
and ip, respectively, and the one-loop effective action S^s reads [33] 



N 



'eflf 



= E 



a,b=l 



3 In sinh 



1 /I 



2 V3 



-/3 + i{aa — cxb) I > + 6 In sinh 



m 



P + i{aa - Oih, 



-8 In cosh ^ ^ (^/? + ^(«a - Oib) 



N 



sin { 2^"" ~ 



In the general case, the one-loop effective action ^eff can be decomposed as 



(3.9) 



eflf ' 



7=1 



(3.10) 



where the first (second) term comes from integrating the square (rectangular) block of the 
fluctuations, which corresponds to the interaction between the "fuzzy spheres" represented 
by matrices of the same (different) size, respectively. One series of the spectrum for Yi 
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gives I Y^/(/ + 1), which coincides, including degeneracy, with the spectrum for the gauge 
field fluctuation A and ghost fields c, c. These arc the unphysical modes, which cancel each 
other exactly. Thus, each term in eq. (3.10) is given by 



4' = E 

a,b=l 

nj 



^ (2^ + 1) Insinh -(- 
1=0 ^ ^ 

+ ^(2/ + 1) In sinh + ^("i'^ - « 

1=1 ^ ^ 

+ ^ 6(2Z + 1) Insinh - ( -P{1 + -) + i{ai'^ - a^' 
1=0 ^ ^ 

n / / 1 \ 

- ^ 4(2Z + 1) In cosh |- i^/3i- + -) + i(ai^) - 4'^)j 

4(2Z + 1) In cosh |i + 1) + i(a(^) - ))) | 



'eff 



ki 

-El- 

fcj kj 

EE 



a=l 6=1 



sin ■ 



(nj+nj)/2-2 



(ni+nj)/2 



+ 



+ 



^ (2Z + 1) Insinh I i + 1) + z{a^^ 

l=\ni-nj\/2-l ^ ^ 

m + 1) In sinh U[^f3l + i{a^P - a^^^ 

-n J 1/2+1 ^ ^ 

-./)/2-l 

^ 6(2Z + 1) In sinh | - l^-f5{l + -) + - a 



(/) _ JJ) 



l=\ni-nj\/2+l 
{ni+nj)/2-l 



(J)- 
b ■ 



=\ni-nj\/2 
(n/+nj)/2- 



Y m + 1) In cosh I i + i) + i(a(^) - 4'^)) 

l=\ni-nj\/2-l 
{nj+nj)/2-l n / / 1 

■ Y 4(2^ + 1) In cosh I - hi- + -) + i(a(^) - -(^^ 

i=|nr-n,|/2+i 



(3.11) 



.(3.12) 



The "In sinh" and "In cosh" terms, which come from bosonic and fermionic fluctuations, 
respectively, yield an attractive potential between the eigenvalues, while the "In sin" term, 
which comes from the Vandermonde determinant, yields a repulsive potential. 

4. Thermodynamic properties of the plane-wave matrix model 

In order to investigate the thermodynamic properties of the plane-wave matrix model, we 
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consider the free energy defined hy F = — Tin Z — const., where the constant is subtracted 
in such a way that hnijv^oo jf? vanishes at T = 0. This convention is motivated from the 
requirement that the free energy F should be of order 1 below the Hagedorn transition [33] . 
Throughout this section we restrict ourselves to the s = 1 case, and use the parameters 
n = ni and k = ki. 

When we fix k in the large limit, the I = 0{N) terms in the effective action (3.11) 
have to contribute in order for the free energy to become of order N'^. This, however, does 
not happen unless /3 = since for /3 7^ the I ^ 1/P terms cancel between bosons and 
fermions. Therefore, the free energy limAr^oo j^r vanishes at any T < 00, which implies 
that the Hagedorn temperature is infinite [33] . Note also that there are only a finite number 
of eigenvalues in the present case, and hence there are no critical behaviors associated with 
the dynamics of the eigenvalues. 

On the other hand, when we fix n in the large N limit, the free energy can be of order 
A^^ at sufficiently high temperature, since the effective action is generically 0{k^) ~ 0{N^). 
When P is large, the attractive potential is insignificant, and the distribute uniformly 
in the interval (— tt, tt]. In this case, we will see that limjv->oo ]fj actually vanishes. As /? 
decreases, the attractive potential becomes more important. At some point the distribution 
of Ua becomes non- uniform and limjv^oo 772- becomes negative. 

This transition, which is interpreted as the Hagedorn transition [33], is associated with 
the spontaneous breakdown of the center symmetry A(t) A(t) + const. 1 at. The Polyakov 
line, which is a useful order parameter for the spontaneous symmetry breaking, is given, 
at the leading order of the perturbation theory, by 




where the expectation value is taken with respect to the effective action (3.11) for the 
gauge field moduli. When the center symmetry is unbroken and the distribution of is 
uniform, we obtain P = in the large N limit. When the center symmetry is spontaneously 
broken, and the distribution of aa is non-uniform, we obtain P ^ . This is analogous to 
the deconfincmcnt phase transition in large N gauge theories [47-49,51,53,56]. 

In what follows we consider the n = 1, 2, 3 cases, for which k = ^ goes to infinity with 
TV. The integration over the gauge field moduli, which play a crucial role in the Hagedorn 
transition, can be done analytically near the transition point. For arbitrary temperature 
we perform Monte Carlo simulation to confirm the analytical results and to obtain explicit 
results in the high temperature regime, which is not accessible analytically. 

4.1 Analytical results near the transition point 

In this section, we investigate the behavior near the Hagedorn transition analytically using 
the eigenvalue density 

k 

a=l 

®This constant is the same for all the classical vacua considered in this paper according to the results of 
appendix C. 
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where 6 £ (— vr, tt]. 

First let us review the analysis for the trivial vacuum given in refs. [33,38]. The 
effective action (3.9) can be written in terms of the eigenvalue density p{6) as 



^eff = 



J J dede'p{e)p{e') |- f; l/„,(/?)e^-(^-^') - in sin Q(e - e'] 



(4.3) 



where we have expanded the "insinh", "in cosh" in (3.11) into Fourier series, and the 
function fm{P) is defined by 

fmiP) = 3e-3'^™ + Ge-^'^™ - 8(-l)™e-2'^"* . (4.4) 

By making a Fourier expansion for the eigenvalue density as 

^ oo 
m=—oo 

the effective action can be rewritten as ^ 

oo ^ 

SeS = V -IPmpIl - fmiP)} • (4.7) 



m 

m=l 



If one of the coefficients of becomes negative, the uniform distribution becomes un- 

stable. This indeed occurs for m = 1 as we decrease /3 from oo. The Hagedorn temperature 
Th = 1/Ph, which can be determined by solving 

l-/i(/3i/) = 0, (4.8) 

is Th 2± 0.0758533. 

Since the I3h = 1/Th turned out to be quite large, we may safely omit the m > 2 
terms in the expansion (4.3). This simplification enables the analytic calculation of the 
thermodynamic quantities by the saddle-point method, which is exact in the large N limit. 
The saddle-point equation for the eigenvalue density reads 

''=^Tem^ /d«V(»'){AW.m(«-.')4co.(l«*-»'))} . (4.9) 
This equation can be solved by using the Ansatz of the Gross- Witten form [58] 



p(e) = J ^ i V ^ - i (1^1 ^ ' (4 10) 

I i\e\>e,,), 



^In arriving at (4.7), we have also expanded the "In sin" term using the formula 

In sin - = — ^ — cosmO + const. (4-6) 



m 
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where ^ci = 2 sin ^ ^J^^ and 



Given the eigenvalue density, the free energy ^ and the Polyakov hne can be calculated as 
(See Appendix D for the derivation.) 

F _ r^. n_ ^^^^^^ 



Ar2 Ar2 \2V2y 2V2 

The above analysis can be extended to the general n case in a straightforward manner. 
The only modifications in the effective action (4.3) are to replace the overall factor iV^ by 
A;^, and to replace the functions fm{P) by 

n—2 n n—1 

f^{P) = ^(21 + l)e-3^"*('+i) + ^{21 + l)e-3^"^' + J^6(2i + l)e-3'3'"('+5) 
1=0 1=1 1=0 

- (-1)^*4(2/ + l)e-^"*(5+^) - (-1)"»4(2Z + l)e-^"*(3+i^) . (4.14) 

' 2 ' 2 

By solving eq. (4.8) in the present case, we obtain the Hagedorn temperature 

' 0.0738901 (n = 2) 
Th I 0.0738526 (n = 3) (4.15) 
^ 0.0738520 (n > 4) . 

Since these values are small, we may neglect the m > 2 terms in (4.3) in the general case 
as well. The expression for the free energy (4.12) gets multiplied by the factor ^ due to 
the modified prefactor in Scs mentioned above, while the expression for the Polyakov line 
(4.13) remains the same. Note that one has to use /i(/3) given by (4.14) for the definition 
of Lo in eq. (4.11). In section 4.3 we will confirm the validity of these analytical results by 
Monte Carlo simulation. 

4.2 High temperature limit 

At high temperature (i.e., small /3), the "Insinh" terms in the effective action (3.11) make 
the eigenvalues attracted to each other against the repulsive force coming from the "In sin" 
terms, and therefore {aa — aj,) is typically of order /?. The free energy can therefore be 
estimated by simply replacing the "Insinh" and "In sin" terms by ln/3, and omitting the 
"In cosh" terms in the effective action (3.11). This gives 

Thus, the leading asymptotic behavior of limTv^oo 777 at high temperature is universal. In 
section 4.3 wc will confirm this asymptotic behavior by Monte Carlo simulation. 
®We have corrected an error in the analytic expression for the free energy given in ref. [38]. 
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4.3 Monte Carlo integration over the gauge field moduli 



In this section we perform Monte Carlo simulation to integrate over the gauge field moduh, 
and obtain exphcit results for the free energy and the Polyakov line at arbitrary tempera- 
ture. 
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0.4 

0.3 

0.2 h 
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0.07 0.072 



trivial 
n=2 



0.078 0.08 



We use the Metropolis algorithm for 
the simulation. At each step we gener- 
ate a trial configuration by replacing one 
of the eigenvalues by a random number 
within the interval (— 7r,7r], and accept 
it with the probability max(l, e"^'^''*'') , 
where AScS is the difference of the effec- 
tive action (3.11) for the trial configura- 
tion from that for the previous configura- 
tion. 

In Fig. 1 we plot the expectation value 
of the Polyakov line near the transition 
point. Our Monte Carlo results are in 
good agreement with the analytic results 
(4.13) including the position of the tran- 
sition point. We have also made a similar 
plot for the n = 3 case, but it turned out to be almost indistinguishable from the results 
for the n = 2 case. In Figs. 2 and 3 we plot the distribution of the eigenvalues. While our 
data agree with the analytic results (4.10) near the transition point, we also start to see a 
small deviation as the temperature increases. 



0.074 0.076 

T 



Figure 1: The Polyakov lino P is plotted against 
the temperature T near the critical point for N = 
2000 in the trivial case, and for k = 2000 in the 
n = 2 case. The dotted lines represent the analytic 
results (4.13). 
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Figure 2: The eigenvalue distribution p{6) 
in the trivial vacuum case is plotted for 
T = 0.075, 0.076, 0.080 and N = 2000. 
The dashed lines represent the analytic re- 
sult (4.10) obtained with the Ansatz of the 
Gross- Witten form. 



Figure 3: The eigenvalue distribution p{9) 
in the n = 2 case is plotted for T = 0.073, 
0.074, 0.080 and k = 2000 . The dashed lines 
represent the analytic result (4.10) obtained 
with the Ansatz of the Gross- Witten form. 
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The free energy can be calculated by 

F = -T{S,s) 



(4.17) 



in the large N limit. In order to obtain the free energy accurately near the critical point, we 
have to make a large N extrapolation. The dominant finite N effects come from the "In sin" 
term in the effective action (3.11) due to its logarithmic singularity when the eigenvalues 
come close to each other. Since the distance between the nearest eigenvalues is of order 
the finite N effects for evaluating the 0(1) quantity 



1 ^ 



iV2 



a^6 



an 



sm ■ 



Ob 



for instance, is O (^7^) • In Fig. 4 we therefore plot the free energy obtained by (4.17) 
at finite N against in the trivial vacuum case for various temperature. Indeed our 
data can be nicely fitted to straight lines, from which we can make a reliable large-A'^ 
extrapolation. The large- limits obtained in this way are plotted in Fig. 5 against the 
temperature T near the critical point. Our results agree nicely to the analytic results (4.12). 
At sufficiently high temperature, we observe that the trivial vacuum gives the smallest free 
energy. As we decrease the temperature, however, it is taken over by the n = 2 vacuum. 
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Figure 4: The free energy F is plotted 

against for N = 500, 1000 and 2000 in 
the trivial vacuum case for various tempera- 
ture. Our data are nicely fitted to straight 
lines. 
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Figure 5: The free energy F obtained by 

the large N extrapolation is plotted against 
the temperature T near the critical point for 
the trivial vacuum as well as for the n = 
2, 3 cases. The curves represent the analytic 
results (4.12). 



In Figs. 6 and 7 we plot the free energy at higher temperature. The deviation from 
the analytic result for the free energy (4.12) becomes pronounced for T > 1 as expected, 
and our data can be nicely fitted to 

F 

where the leading asymptotic behavior is determined analytically in section 4.2. Although 
the leading term is the same at A'^ = 00, the coefficient ci of the linear term turns out to 



-8riog(T)-cir-c2 



(4.18) 
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be larger for the trivial vacuum case. Hence we conclude that the trivial vacuum gives the 
smallest free energy for T > 0.077. 




2 4 6 8 10 

T 



Figure 6: The free energy is plotted 
against T for N = 500 in the trivial vacuum 
case. The solid line represents the result for 
fitting the data for 1 < T < 10 to (4.18), 
where a = 15.852(1) and C2 = 0.038(3). The 
dotted line represents the analytic result for 
the free energy (4.12). 




-350 I ' ' ' ' 1 

2 4 6 8 10 

T 

Figure 7: The free energy is plotted 
against T for n = 2 and k = 500. The 
solid line represents the result for fitting 
the data for 1 < T < 10 to (4.18), where 
ci = 11.948(5) and C2 = 0.16(1). The dotted 
line represents the analytic result for the free 
energy (4.12). 



5. Summary and discussions 

We have studied the thermodynamic properties of the plane-wave matrix model in general 
classical vacua. The Hagedorn transition occurs due to the dynamics of the gauge field 
moduli. Extending the previous works on the trivial vacuum, we obtained analytical results 
for the thermodynamic quantities near the transition point, which revealed the existence 
of the "fuzzy-sphere" phase. We also performed Monte Carlo simulation to integrate over 
the gauge field moduli. This method can be used for arbitrary temperature. Usually free 
energy is difficult to calculate by Monte Carlo simulation, but in the present case we can 
obtain it very accurately from the expectation value of the effective action for the gauge 
field moduli, thanks to large N. We observe that the trivial vacuum gives the smallest free 
energy in the high temperature regime. 

The plane- wave matrix model is closely related to A/" = 4 SYM theory on i? x S*^ [59] , 
which implies that the matrix model at finite temperature is related to the SYM theory 
on X . We therefore expect that our results should have implications on the phase 
structure of the SU{oc) SYM in four dimensions. In particular, the fuzzy sphere phases 
we found in this paper may correspond to the plasma-ball phases in the SU{oo) SYM on 
X S*^, which are interpreted as localized black holes on the gravity side via the AdS/CFT 
duality [50]. Another interesting future direction is to clarify the relation between our result 
and the entropy bound discussed recently [60]. 
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We hope that our work provides a clue to the phase structure of the AdS black hole 
and to a deeper understanding of the AdS/CFT correspondence at finite temperature. 
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Appendix 

A. Gauge fixing 

In order to arrive at our main result eq. (3.8) in section 3, we have first integrated over 
the fluctuations fixing the gauge field moduli, and then the integration over the gauge field 
moduli has been written in terms of the moduli parameters ai^ . Since the classical vacua 
(2.4) breaks the U{N) gauge symmetry down to 11/=! U{ki), we have to make gauge fixing 
at each step of the above procedure. First, when we integrate over the fluctuations, we 
have to fix the gauge along the gauge orbit in the direction of the coset space rr^^^^'TnT^- 
This is necessary since otherwise there will be zero modes corresponding to the broken 
symmetry, and we cannot integrate over the fluctuations. Next, when we integrate over 
the gauge field moduli, we have to fix the gauge corresponding to the remaining symmetry 
YYi^i U{ki) in order to reduce the integration to that over the moduli parameters a^P . This 
is necessary for studying the large N limit analytically by the eigenvalue density method 
and for simplifying Monte Carlo calculations drastically. In what follows we explain the 
two steps of gauge fixing procedures separately. 

A.l Integrating the fluctuations 

Let us explain the first gauge fixing, which is necessary for integrating the fluctuations. 
We use the background field gauge 



where r{t) is an arbitrary function, and and A^(t) for = t, 1, 2, 3 are defined by 



D^A^{t) = r{t) , 



(A.l) 



Dt = dt- i[A{t), ■ ] , At{t) = A{t) 
Di = -i[Bi, ■], A{t) = X,{t) 



for i = 1,2,3. 



(A.2) 
(A.3) 



The corresponding Faddeev-Popov (FP) determinant is defined by 




(A.4) 
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where g{t) is an element of the coset groTip Yl'^^^jlkj) ^m(^) represents the gauge 
transformed field. The FP determinant can be represented as 

A^^l = J [dc] [dc] exp j dt tr [DfC ■ DfC - [Bi, c] [Xi, c])^ , (A.5) 

where c, c are the ghost fields satisfying the periodic boundary condition. Inserting the 
identity (A. 4) in the partition function, and integrating over r{t) with the gaussian weight 
exp (— 5 / dtr'^{t)), we obtain the gauge fixing term and the ghost term 

5g.f. = J dtti 1^ (D^A,,)^^ , 5ghost = J dtti {Dtc-Dtc-[Bi,c][Xi,c]} , (A.6) 

which should be added to the original action (2.1). 
A. 2 Integrating the gauge field moduli 

Let us explain the second gauge fixing, which is necessary for integrating the gauge field 
moduli. We use the static diagonal gauge, which is analogous to the one adopted for the 
trivial vacuum [33]. 

First we impose the static gauge condition given by 

dtA^^\t) = 0. (A.7) 

The corresponding FP determinant can be defined by 



where g^^\t) is an element of the unbroken gauge group 11/=! U{ki). The FP determinant 
App is given by 

s 

A.^^l = \[dei'(dtD^P) , (A.9) 



1=1 



where we have defined D^^ = dt — i[A^^\ ■ ], and the symbol det' implies that we have 
omitted the zero modes. The static gauge (A.7) does not fix the gauge completely, and we 
still have the global 11/=! U{ki) symmetry as the residual gauge symmetry, which we fix 
by further imposing the constant matrices A^^^ to be diagonal (3.2). 

As is well known in matrix models, this gauge fixing yields the Vandermonde deter- 
minant, which is derived as follows. The path integral measure for the gauge field moduli 
[(i^(i)] around the static diagonal configuration (3.2) is rewritten as 



s I . ki [I) k, ki \ 

n o ■ n ^ ■ n w - -i" ■ ■ n n<:.. . (^-^ 



, kjl -L-L 27r 

7=1 Y a=l a^b m^O a,b 



where g^^'^ is an clement of n/=i and dAm {tti ^ 0) is the integration measure for 

the non-zero Fourier modes. 
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Let us evaluate the determinant (A. 9) explicitly. In terms of Fourier modes ^, the FP 
determinant (A. 9) reads 



A 



(2) 
FP 



n 

7=1 



Yl n (27rimr)|27rimr + i{ai^^ - ^)r} 



a, 6=1 mj^O 



(A.12) 



Therefore the identity (A. 8) is rewritten as 

kr 

a,b=l m^O 



1=1 



mjtO 1=1 



kj oo 



dg"'\t)YlYl 

ajtb m=l 



27rm 



(A.13) 



In the second line, we have omitted the delta function for the non-zero Fourier modes, 
which are integrated out by using the measure in (A. 10). Inserting the identity (A.13) in 
the partition function, the measure (A.IO) for the gauge field moduli [(i^(t)] can be written 
as 



[dA{t)] = (2r)^^=i'=? • Yl (27rimr)S^=i'=? 

s ( -I ''I J (-f) 



1=1 



ki\ -L-L AttT 

a=l a^b 



sm 



(A.14) 
(A.15) 



where dg^^\t) = dg^^\t) ■ dg'^^^ represents the measure for the gauge function of the 
U{kj) symmetry. The "sin" term in (A.15) appeared from the formula used to sum 
over the Fourier modes 



n > 



n=l 



n 



X" \ SmTTX 
2 



TTX 



The result for the trivial vacuum can be obtained formally by setting s = 1 , fci = A'^ 
and ni = 1 , which yields 



(2Tr.n(2«».T)«'^<i.(o.in^^n 

m^O a=l a^b 



sm ■ 



aa - ab 



(A.16) 



where g{t) represents the gauge function of the U{N) symmetry. This agrees with the 
measure obtained in ref. [33] . 



^The fields M{t) and Q{t) obeying the periodic and anti-periodic boundary conditions, are expanded as 



(2m-l) , 



m= — oo 



(A.ll) 



m= — oo 



respectively. 
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B. Evaluation of the determinant 



Integrating over the fluctuation Y , A and ghost fields using the quadratic action (3.6), 
we obtain the determinants 



detppt + A) , 
c,c: detB(A^ + A2), (B.l) 

where A represents an eigenvalue of the mass operators (3.7), and the lower suflfices of det 
specify the boundary condition (periodic for "B" and anti-periodic for "F") for the fields 
on which the operators act. 
Using the formulae 

sinhvra; -tt / \ ■, /vrx\ 



n=l ^ ' n=l 

we obtain the determinants as 



det [-d] + >?) =\{\{ {2TrimT + i{aa - ab)T + A} 

m a,b 

= {2Tf' ■ n {2'KimTf" -JlsinhWipX + iiaa-a,))] (B.2) 



12 

m^cO a,b 

det F ( A + A) = n n {^^(2™ - 1)^ + ^("a - cib)T + A} 

m a,b 

= n (^^(2m - l)r)^' • II cosh I ^ (/3A + i(a„ - «&)) | . (B.3) 

m a,b ^ 

The coefficients in front of "sinh" and "cosh" have to be taken into account in comparing 
the free energy for different vacua. 

C. The overall factor of the partition function 

In this section we obtain the overall factor C in (3.8). There are three types of contributions 
to C, which affect the calculation of free energy at the leading order in N; namely the 
Gaussian integration, the determinants, and the gauge volume obtained from gauge fixing. 
The factor obtained from each contribution is given in table 3. The symbol Vol(C/(A^))iocai 
represents the gauge volume for the local U {N) symmetry, and we have introduced the 
functions 

/(T) ^2T H {2mmT) , ^(r) = n{vri(2m - l)r} . 

m^O m 

The overall factor of the partition function C turns out to the same for all the classical 
vacua, and it is given by 



9Ar2 



C ^ n 2^ • /(T)-^^' • giTf"" ■ Vol(f/(iV))io,„ . (C.l) 
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Gaussian 


determinant 


gauge volume 


Jdgit) 






Uif=iVoi{(7(fez));,_, 








n/ilV0l(f/(fc/))local 


r{t)-mt. 








y(i)-int. 








A{t)-mt. 




j-(y)-(iv2-Ef=ifc?) 




<li{t)-mt. 








c{t), c(t)-int. 








[<IA] 








C 






Vol(f/(iV))iocal 



Table 3: The list of contributions to the overall factor of the partition function for general classical vacua. 



D. Evaluating free energy with the eigenvalue density 



In this section we derive the free energy (4.12) for the eigenvalue density (4.10). Omitting 
the m > 2 terms in (4.3), the effective action is given as 



-TV' 



6 fO 

'^1^1 dede'p{e)p{e') h{f3) cos {e -e') 



-N^V.V. 



ded9'p{0)p{e') In 



sin[-{9-e') 



1 r 



dein 



sm ■ 



(D.l) 



where the symbol V.V. represents the principal value. For eigenvalue density (4.10) one 
can easily obtain 



J-0C 

r0ci 
V.V. / 



sin 6* = , 



dOpiO) cos 6* = 1 



dein 



9 



sm ■ 



d9d9'p{9)p{9') In 



sin ■ 



■ln2 + - In 

2 2 4 



■In 2 (D.2) 
(D.3) 



where, in deriving the formula in the second line, we have used the fact that the distribution 
(4.10) satisfies the saddle-point equation (4.9). Using these formulae, we obtain 



-^«/i(/?)(l-4j +21^2-4 



, 1 u 1 
2--]+-ln--- 



(D.4) 
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